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ABSTRACT 

We find the cross-over behavior for the spin-spin correlation function for 
the 2D Ising and 3-states Potts model with random bonds at the critical point. 
The procedure employed is the renormalisation approach of the perturbation 
series around the conformal field theories representing the pure models. We 
obtain a crossover in the amplitude for the correlation function for the Ising 
model which doesn't change the critical exponent, and a shift in the crit- 
ical exponent produced by randomness in the case of the Potts model. A 
comparison with numerical data is discussed briefly. 
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1 Introduction 



In the studies of critical phenomena in real physical systems, impurities and inhomo- 
geneities are always present. Many theoretical models have been proposed for the study 
of random models. In most of these cases the replica trick for quenched systems is em- 
ployed. This corresponds to computing the averaged free energy by taking n copies of the 
system and going to the limit n — 0. For models with random bonds, the main problem 
is to determine if the randomness leaves unchanged the critical properties of the pure sys- 
tem or if the singularities of the thermodynamical functions are eliminated. First results, 
which suggest an intermediate situation, have been obtained by Harris and Lubensky [|l|, 
Grinstein and Luther and Khmelnitskii using the standard 0^ theory. Other cases, 
like long-range correlated quenched defects have also been considered, (see for example 
IQ, 0.) A first step in the understanding of the relevance of randomness was given by the 
Harris criterion 0: the randomness is relevant (irrelevant) if the specific heat exponent 
of the pure model is positive (negative). Two dimensional systems are particularly inter- 
esting because of the rich structure of conformal invariance in this dimension. Assuming 
that a random model has a critical point with second order phase transition, the main 
interesting problem is to determine which conformal field theory represents this model 
at the infrared fixed point. Let us also mention that an exact result has been obtained 
by McCoy and Wu who considered a two-dimensional Ising model where only vertical 
bonds on a square lattice were allowed to acquire the same random value. They found 
that the logarithmic singularity of the specific heat disappeared completely. 

The models that we will study in this paper are the two dimensional Ising and Potts 
models with random bonds. For the case of Ising model, Harris criterion doesn't provide a 
qualitative answer of the relevance of the randomness (because the specific heat exponent 
for the 2 — D Ising model is 0). First results were obtained by Dotsenko and Dotsenko 
1^ who showed that near the critical point, this model can be represented by an n = 
Gross- Neveu model 0. With this technique, they found that the specific heat singularity 
get smoothed as Zn(Zn(-pj)) where \t\ is the reduced temperature. Calculation of spin-spin 
correlation function by this technique which involves non-local fermionic representation of 
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< o"(T >, was later questioned by Shalaev and Shankar [|T^, 11] who gave arguments that 
the asymptotic behavior of this correlation function is unchanged by the randomness (see 



also [T^.) Some arguments on why computations using non-local fermionic representation 
fail to give the correct result for the spin-spin correlation function in the random case will 
be given at the end of the paper. More recently, this situation was questioned by Ziegler 
|T3| . This author claims that non-perturbative effects introduce an intermediate phase 
around the critical point of the pure model. On the other hand, numerical simulations of 



the Ising model [p!5| , seem to confirm the theoretical predictions of the specific heat 
and spin-spin correlation function asymptotic behavior. For the Potts model, Harris 
criterion predicts that randomness is relevant and changes the critical behavior. Using 
conformal field theory techniques, Ludwig also perturbatively computed a shift in 
the critical exponent of the energy operator in the case of the random Potts model (this 
critical exponent for the Ising model is also unchanged). 



The paper, which complete the work presented in is organized as follows. In 
section 2 we explain how the replica trick method can be used to deal with the parti- 
tion function of the 2D Ising and Potts models with random bonds. Near the critical 
point, these models are represented by perturbed conformal field theories. We give a 
brief summary of the Coulomb Gas representation of conformal field theories needed to 
pertubatively compute correlation functions. We recall also a kind of "e" regularization 
which consists of shifting the central charge of the pure model. Then, in section 3, we 
sketch the renormalisation group equations for the coupling constant and correlation 
function. Specializing then to the particular case of the Ising model, we find a cross-over 
in the amplitude at intermediate distances produced by randomness. In the case of the 
Potts model, we show that the spin exponent is modified in a power series of e due to 
the randomness. Finally, section 4 contains a summary of our results and a discussion. 
Technical details of the renormalisation group calculation and the calculation of integrals 
are put into Appendices. 
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2 The model 



Our starting point is the 2-D Ising-Potts Hamiltonian at the critical temperature: 

•^o = E i(^J^)K,s, (2.1) 

<i,3> 

where Jo is the couphng between spins, Si are the spin variables and < i, j > means that 
the sum is restricted to nearest neighbor spins. The unperturbed partition function is: 

Here, Tr^. means a summation over all spin configurations. The addition of a posi- 
tion dependent random coupling constant produces the following change in the partition 
function: Jq ^ Jq + AJj It is well known that at the critical point the unperturbed 
models are represented in the continuum by unitary minimal conformal field theories in 
which the relevant operators are the spin operator a corresponding to Si and the energy 
operator e corresponding to 5s^,sj ■ In the continuum the term A Jj ^ can be written as a 
position dependent mass term m{z), which gives the total partition function in the form : 

In order to compute the free energy, we use the replica method. Taking the partition 
function of n identical copies of the system and analytically continuing to the limit n — > 
gives the quenched free energy 



-I3F = ln(Z) = lim 



- 1 



n-»0 n 

where: 

n . n 

- Yl ^o,a - / rn{z) E Sa{z)(fz 
= n Tra,s,e «=i »=i (2.3) 

o=l 

the average of is made with a Gaussian distribution for mlz): 

Z^^ [l[dm{z)Z^e-^^"'^'^-"'°^" 

•' z 

which gives: 

n . n . n 

- E 'S'o,a + 9o J Y Sa{z)eb{z)(Pz -mo J Y ^a{z)(fz 

a=l 



The average with a more comphcated distribution gives higher order cumulants, produc- 
ing in ( ^.41 ) terms which can be shown to be irrelevant by power counting Also, 

n 

the terms in / ^ ea{z)eb{z)(f z containing the same replica label produce irrelevant 

a,b=l 

operators and can be omitted. We obtain then the following partition function: 

n „ „ n 

t[Tra,s,e ^^=1 (2.5) 

a=l 

In the limit mo — > 0, this model corresponds to a conformal field theory perturbed by the 
term that is quadratic in the e operator. Then the "evolution" of the coupling constants 
Qq and mo under a renormalisation group (R.G.) transformation can be analyzed as well 
as the behavior of the correlation functions. In the calculation of correlation functions 
< 0{0)0{R) >, where O is some local operator, we will proceed perturbatively: 

< 0{0)0{R) >=< 0{0)0{R) >o + < SiO{0)O{R) >o +i < S]O{0)O{R) >o +... 

where <>o means the expectation value taken with respect to 5*0 and 

Si = J Hj{z)d^z = gJ Y.^-(^)^b{z)d^z (2.6) 

The operator O is then renormalised as 

O ^ 0(1 + A^g^ + A^gl + ^3^70' + ■■■) = Z^O 

The integrals of correlation functions involved in the calculation can be performed by 
analytic continuation with the Coulomb-gas representation of a conformal field theory ||T8| 
where the central charge is c = |+e'. The e' term corresponds to a short distance regulator 
for the integrals. In addition, we also used an infrared (I.R.) cut-off r. The result is then 
expressed as an e' series with coefficients depending on r. The limit e' ^ corresponds 
to the pure Ising model at the critical point while the Potts model is obtained for some 
finite value of e'. We recall here some notations of the Coulomb-gas representation for 



the vertex operators |]T8|]. The central charge c will be characterized in the following by 
the parameter = 1^-^ = | -|- e with 



1 - 24a^ ; a± = ao± J al + 1 (2.7) 



Note that for the pure 2D Ising model = | and c = | while for the 3-state Potts 
model = |, c = I and e = — ^. The vertex operators are defined by 

Vn^{x) = e--<^(-) (2.8) 

where is a free scalar field and where the a„m are given by 

a„m = ^ (1 - + ^ (1 - m)a+ (2.9) 

The conformal dimension of an operator Vnm{,x) is A„m = —oinmOinm with 

amn = 2ao - a„m = ^ + ^)"- + ^ + "^)"+ (2-10) 

The spin field a can be represented by the vertex operator \4>,p-i whereas Vi,2 corresponds 
to the energy operator e. In the same way, we associate e*°+'^'^^') to the screening charge 
operator l + (x). Note that in the Ising case the a operator can also be represented by 
the V21 operator (since both operators coincide in the limit e — 0). So, we can represent 
our spin operator by Vk,k-i where k = 7^^- We have A = 2 for V21 and A = | for 



3 Renormalisation Group Equations 

In this section, we will deal with the computation of correlation functions of operators e 
and a. To compute them, one needs to determine the effect of the random coupling on 
the operators e and a and compute the renormalised operators e' and a'. This means 
that we want to compute the functions and Z^- such that 

e' = Z,e and a' = Z^o (3.1) 

A convenient way to define and Z„ is to consider the more general action 

n „ n . n . n 

Y.^o,a-9o J £a{z)eb{z)d'^z + mo I ^a{z)d'^z - ho I Y^a{z)d^z (3.2) 

a=l a,h=l a=l a=l 

This merely corresponds to the action used in ( |2.4| ) with an additional coupling of the a 
field. Then, with the help of the operator algebra (O.A.) coming from contractions 
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between e and a operators, we will compute the effect of the fifo / ^a[z)eb{z)d 

a,6=l 



term (i.e. the random part of the model) on the coupling terms mo / ''^ea{z)^z and 

n 

/lo / X/ ^a{z)'^z. More precisely, we will compute 



0=1 



a=l 



ffo / £a(^)£b(^)c?^^ I rxiQ f f2^a{z)d^z \ ~ m / J2^a{z)d'^z (3.3) 

j \ V a,&=l J a=l I •' a=l 

and 

E h?oy E ^a('^)^f('^)^''^ Poy E^a('^)^''^ C^hj J2aa{z)d^Z (3.4) 
« y \ a,&=l / a=l J a=l 

m and h being the renormalised coupling constants. Obviously, this computation will 
be perturbatively made only up to some finite power in qq. In fact, the first step of 
the computation will be to determine the renormalised g constant on which and Z„ 
depend. 

3.1 Renormalisation of the coupling constant g. 

The renormalisation of the coupling constant g will be determined directly by a pertur- 
bative computation, g is also given by the O.A. producing 

9oJj2 £a{z)£b{z)d^z + ^IgojYl £a{z)£b{z)d'^z\ + (3.5) 

a,b=l \ a,b=l j 

+ M^oy E £a{z)eb{z)d^z\ H = ^y" ^ ea{z)eb{z)d^z 

\ a, 6=1 / a, 6=1 

with g = g^s^ ^290 + ^3^0 + " " " where A2 comes from 

ea{z)eb{z)d^z j £a{z)£b{z)d^z ^ j £a{z)£b{z)d^z + • • • (3.6) 

a,6=l a,b=l a,b=l 

and from 



.2^ 



g y E £a{z)£b{z)d^z j Y £a{z)eb{z)d^z j Y £a{z)sb{z)d^ 

a,b=l a,b=l a,b=l 

„ n 

^A^j Y ea{z)eb{z)d^z + --- (3.7) 



a,6=l 



Here, we will only perform the computations up to the third order in g^. The technical 
details of the computations of A2 and are given in the appendix A with the following 
result for g{r) : 

9{r) = r-^' {90 - gl^^n - 2)^ + gl^^^n " 2)^ (l + ^^^)) (3-8) 

Note that we multiply the result by r~^^ in order to obtain a dimensionless coupling 
constant g{r). 

In the computation of correlation functions, we will need the /9-function associated to 
the renormalisation group equations of g{r). It can be derived directly from the previous 
expression of g{r) with the result : 

m = ^ = -3e^(r) + 47r(n - 2)g\r) - l&^\n - 2)g\r) + 0{g\r)) (3.9) 

Finally, taking the limit n ^ 0, we obtain for the /9-function up to the third order : 

= -3e5 - 8V + 327rV (3-10) 

We can then immediately note that, in the limit e ^ {i.e. the random Ising model), 
j3{g) has an infrared fixed point aX g = {). In the case where e < {i.e. the random Potts 
model) the infrared fixed point is located aX gc = + ^ + 0{e^). 



3.2 Renormalisation of a and e 



In order to be able to compute the correlation functions of a and £, the second step is 

to determine the effect of the renormalisation on these operators. One needs to compute 
the multiplicative functions Z^j and Z^. This will be made by computing the renor- 

n ,. n n 

malised coupling terms m J '^ea{z)(fz — {rrioZ^) / ^ ea{z)(fz and hj^2 <^a{z)d^z — 

a=l a=l a=l 

/n 

aa{z)<f'z. A direct computation of both m and h will provide us with the 
a=l 

functions Z„ and Z^. As for the computation oi g — Zggo, we will compute in perturba- 
tion : 

"^0 j Y,^o.{z)(fz +gomo j ea{z)eb{z)d^z j ^ea{z)d^z+ (3.11) 

a=l a, 6=1 a=l 
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,9l 



/ ^ ea{z)e^{z) / ^ ea{z)d^z H = m / ^ £:a(^)c?^ 

^•^ a,fe=l / a=l •' a=\ 

and m = mo(l + -Bififo + -B25'o + ' ' ') with defined by 

a, 6=1 a=l a=l 

and B2 by 

I ea{z)e,{z)\ Jj2^aiz)dh = B2jY^aiz)d' 

y a, 6=1 J a=l a=l 

The details of the computation are presented in the appendix B, with the result 
r-^+2^m(r) = mo(l - 47r(n - l)^7o^ + 47r2(n - l)^7o^(l + ^^)) 

Here again, we multiply m(r) by r^^+i*^ in order to obtain a dimensionless coupling 
constant. We also give the R.G. equation for Z^, which we will need later. 

dln{Z,{r)) 



z 



dln{r) 



A7r{n - l)g - 8TT^{n - l)g^ 



(3.12) 



Similarly, for the coupling constant Hq, we compute up to the third order 
J Y ^a{z)eb{z) J Y.^a{z)d^z = C, j Y(^a{z)dh 

y a, 6=1 J a=l a=l 

We give here directly the result (see eq.( |C.7| ) in appendix C) : 



r-f-''(^)h{r) = /io(l + (n-l)^oV 



-6e 



-12{n-l){n-2)gl^ 



^^P(-|)P(-|) 



(3.13) 



9e 



l + 9(2-A) 



r^(-i)r^(|) 
r^(-|)r^(-|; 



The multiplicative term r"^"'^^'^'' in front of h{r) is introduced in order to make this 
parameter dimensionless. Here, a(e) is a function of e depending on which representation 
of the spin field we are taking in the Coulomb gas picture (see section 2). Its explicit 
form will be irrelevant in the following. The corresponding R.G. equation for Z^- will be 
given by 



dln{Z„{r)) 
dln{r) 



-3{n-l)g^{r)n^e 
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3' 



r'(-i)r'(i 



+4(n- l)(n-2)7r^^^(r) 



(3.14) 



4 Correlation Functions 



We now have all the ingredients needed in order to compute the correlation functions. 
They will be calculated with the help of the R.G. equations, for the theory with mo, 
Hq ^ 0. From the R.G. equations, we have : 



< e{0)e{sR) ^^^^lifL^s''''^ < e{0)e{R) 

Z^[r, g[r)) 

This can be written as : 



< £(0)e(si?) >r,g(r)= 6 5"^^^ < £(0)e(i?) >r,g(^sr) (4.1) 

where we used the notation : 

dlnZg , , / . „\ 

-TT^ = leig) (4.2) 
alnr 

and (^(s) = g{sr)]gQ = g{r). We assume now r to be a lattice cut-off scale. In a similar 
way for < a{0)a{R) > the R.G. equation is : 



2 f 2^dg 
J 13(g) ^ 

< a(0)cr(si?) >T-,g(r)= e ''O S " < cr(0)cr(i?) >r,g{sr) (4.3) 



with 

dlnZ„ , , , , 

In equations ([4.1| )- (^4.3| ), R is an arbitrary scale which can be fixed to one lattice spacing 
r of a true statistical model. The dependence of < cr(0)(T(r) >r,g(s) on s will then be neg- 
ligible, assuming that there are not interactions on distances smaller than r. Therefore, 
it reduces to a constant. Then, s will measure the number of lattice spacings between 
two spins in < a{0)cx{sR) >. In the following, when we treat separately the Ising model 
and the 3-states Potts model, we adopt the choice r = 1. 

4.1 The Ising model 

The Ising model corresponds to the case e — and so the (3 function is : 

(3{g) = -8ng^ + 327rV (4-5) 
10 



Therefore, we can see that the I.R. fixed point is located at g = 0. Also we have, by 
eqs.(g3g), for n = 0, e = and definitions Q, 



je{g) = -47r(7 + 87rV (4.6) 
7.(^7) = 8ng^ (4.7) 
The integral for the e correlation function, eq.( [4.1| ), gives : 

9{s) 9{s) g{s) 

50 ^^^^ 90 ^ ^ 90 ^ 

= 2n{gis)-go) + ln(^^^ (4.8) 
Now, we need to compute g{s). The integration of equation /3{g) = -^^^ gives : 

a{s) , sr 

rin r 

dlnr 



-87r^2 _^ 327rV 



90 



with the following solution up to the second order : 

( \ 30 , 47r^o/n(l + 87r^o/n(3)) 3 

= l + 8.golnis)^^ + 1 + SngoHs) ^ + ^^^^^ ^'-'^ 

So, < > correlation function is given by : 

< siO)sis) ^ (1 - 2n{go - g{s))) s'^^^ (4.10) 

9o 



For the cr correlation function, the computation is similar. In fact, in that case, we have 
•jaig) = SiT^g'^. Thus, keeping only the first order of the /^-function {i.e. [3{g) = —Sng^), 
we obtain : 

9{s) , , 9(s) 
lA9) 



2 / ^dg = -2vr2 / gdg = -rr' (gis)' - g^) + 0{gl) (4.11) 

90 ^^^'^ 90 

The < aa > correlation function is then found to be given by : 

<a(0)a(s) > ~ (1 + TT^ - SW')) s-^^- (4.12) 

11 



The calculation of the term in the j3 function and the g"^ term in the renormalisation 



of e was already done in |TB[, extending the one loop result of [^. We recovered these 
higher order corrections using a different technique, which allowed us to calculate also 
the modified correlation function of the spin operators. 



4.2 The Potts model 



We now consider the 3-state Potts model. With the convention that we use this case 

2 



corresponds io t = (3{g) is given in eq. (|3.10| ), and 

lei.9) = 



-Ang + Sn'^g'^ 



Svr^e l + -(2-A) 



r^(-|)r^(|; 



(4.13) 
(4.14) 



. 3^- '•^n(-l)P(-i)^ 
At long distances, the integrals in eqs. ([4.1| ), ([4.3|) will be dominated by the region g ~ gc, 
with fi'c = — If: + + 0{e^). This is different from the Ising model, because here '~ie{9c) 
and 7o-(fi'c) have finite values for g = gc. Thus, 



9is) 



dg ~ le{gc)ln{s) 



and 



90 ■ ' go 

The correlation functions can then be deduced directly : 



ia{g) 



dg ^ '^e{gc)ln{s) 



< e{Q)e{s) >g, ~ s-(2Ae-27.(9c)) < a(0)(T(s) >g„~ s 



(4.15) 



(4.16) 



So, we can see that a direct consequence of the new IR fixed point is a modification of 
the critical exponents and A^.. A straightforward computation will give these new 
exponents : 



2„2 

c 



2A; = 2A, - 27,((7c) = 2A, + 87r^7c - 16 A 

9 
4 



2A, - 3e + -e^ + O(e^) 



(4.17) 



and 



2 a; = 2 A, - 27, ((7e 



2A.-3(2-A)^A^,3^0(e^) 



- iQ^'gl 



(4.18) 
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Here A = | and the final result for the new critical exponent is : 

27 r2('_2\p2('l\ 

^K-2A,--^A^A^^,^0i.^) (4.19) 

The value of the critical exponent for the energy operator was already computed, up to 
the second order, by Ludwig |TB|. For the spin operator, deviation of the critical exponent 
from the pure case appears only at the third order, which we compute. This deviation is 
in fact very small. While in case of 3-state Potts model without disorder 2A^ = ^, we 
obtain the new critical exponent 2A; = ^ + 0,00264 = 0,26931. Thus, the deviation 
corresponds to a modification of 1%. 



5 Conclusions 

In case of Ising model spin-spin function the calculation of up to third order of the 
renormalisation group was needed to find the deviation from the perfect model case, in the 



form of the cross-over in the amplitude. This completes the observations of ||10|, [IT], p!2| , 
based on absence of renormalisation of this function in the first order, that asymptotically 
the spin-spin function has the same exponent as in case of the perfect lattice model. 

Earlier calculation of the spin-spin function in the second reference of Q, in which 
a non-local fermionic representation had been used, gave a different result, which was 
not confirmed neither by further theoretical calculations, based on local renormalisation 



directly of spin operators, nor by the numerical simulations in |T^. In calculation 
of 1^ a double summation series had been involved, both sums being log divergent. 
One sum was due to non-local representation of the spin-spin function, another was due 
to randomness- caused interactions. In the calculations, a particular ordering of the two 
summations had been employed : renormalizing the terms of the first series and summing 
up next. We may guess, lacking a better established reason, that such a treatment might 
not be justified. 

Recently the numerical simulations of the random Ising model has been performed 
which measure directly the deviation of < aa > from the pure Ising model at the critical 



13 



point |jT5[|. These measurements were made for disorder such that SngQ ^ 0.3. Deviations 
predicted by our computations are very small. They correspond to 0.1%. The deviations 
obtained in numerical simulations are around ten times larger, and they are of opposite 
sign, i.e. correspond to an extra decrease of the spin-spin function with distance r. 
In it has been checked that this decrease corresponds, within the accuracy of the 
measurements, to a factor function of the ratio r/L, F{r/L), r being the distance between 
the spins and L is the lattice size. So they correspond to finite size effects, being different 
for perfect and random models. We would suggest, on the bases of our calculation of the 
r dependence of the spin-spin function on an infinite lattice, that numerical deviations 
will continue to be plotted by the same curve F{r/L), if one measures < aa > for 
different lattice sizes as it has been done in [|T3|, until the accuracy reaches the value of 
the r-deviation which we calculated here. Only then the curves for different L will split. 

In case of Potts model, the deviation of < cro" > appear also only in the third order like 



for the Ising model, this time in the exponent of the asymptotic function, eq.( |4.19| ). This 



might be easier to check with numerical simulations [^. We refered to the 3-state Potts 
model in Sec. 4. 2 because this particular case might be easier to realize numerically, and 
also because we expect that e-expansion should be well defined in this case. Extending 
our result to the 4-state Potts model might be questionable, as this is actually the limiting 
case in the range of minimal conformal theories which are parametrized by e. 
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A Renormalisation of g 



In this appendix, we shall describe the calculations for the renormalisation of g up to the 
second order. 

A.l First order 

The first order correction to g comes from the contraction of two (ee) terms (plus some 
combinatorial factor corresponding to all of the possible contractions) : 

2 

\x-y\<r '^^'^ 

When b — c ^ a,d this contractions gives : 



Here, A2{r, e) is the result of the integral produced by the contraction of two e operators : 
A,{r, 6) = 2(n - 2) J < e{x)e{y) >o d'y = 47r(n -2) J l-^^) 

\y-x\<r y<r 

= -47r(n - 2) j (A.l) 

< e{x)e{y) >o being the standard correlation function of the unperturbed conformal field 
theory. 

A. 2 Second order 

We will now compute the second order corrections to g, produced by the contractions of 
three {ee) terms. Here, there will be two contributions coming from : 

^ / J2^'^i^)^b{x)J2^c{y)Ed{y)J2^e{z)ef{z)d'^xd'^yd'^z (A.2) 

■ \x-y\,\x-z\<r "^'^ 

The first one corresponds to the case when b = c; d = e; b, d ^ a, f . We will denote this 
term by : 

A3^i{r,e)gl / J2 ^a{x)ef{x)d'^x 
15 



where >l3,i(r, e) corresponds to the following : 

^,i(r, e) - 4(n - 2)(n - 3) J < e{x)e{y) >o< e{y)e{z) >o (fz£y (A.3) 

\y—x\,\z—x\<r 

Now, replacing < e{x)e{y) >o by \x — y\~^~^'^ and performing a trivial change of variable, 



this gives : 



^3,i(r,e) = 87r(n-2)(n-3) j j \z\-'-'^\z - l\-'-'^d'z 

y<r 

= 167r^(n-2)(n-3)(^) (A.4) 
The second term is produced when a — c — e^f]h — d^f. This will be denoted by : 



with: 



^3,2(r, e) = 4(n - 2) y < e{x)e{y)e{z)e{Qo) >o< >o d^yd^z (A.5) 

|2/-a;|,|z-a;|<r- 

Here, the < e{x)e{y)e{z)e{oo) >o term corresponds to the result of eee e obtained by 
projecting eee over 5(00). In order to compute A3^2{i^,^), we need to use the Coulomb 
gas representation for this 4— points correlation function : 

^3,2(r,e) = 4(n-2)Arx 

I < V-^{x)Vu{y)Vuiz)V+iu)Vuioo) >\y- z\-^-^'d^ydhd^u (A.6) 

\y-x\,\z-x\<r 

N is a, normalization constant which can be fixed by the OA relation in the following 
way. We have : 

< s(0)s(R)s(x)e(y) >^ N J < Vj2{0)Vu(R)Vu(x)V+(u)Vi2(y) > d\ 



In the limit i? — > 



( / is the identity operator) and 

\x - y\ 



< e{Q)e{R)e{x)e{y) < e{x)e{y) >= 
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For the r.h.s. we obtain 

J < V2aMVuix)V^,iy)VM > d'u 

which produces 

i?4Al2|a;_y|4Ai2 y I I I I 

Comparing both resuhs we obtain : 

N = U wr^^'^-w-ir^^-d^-' = -\^-^^ (A.7) 



v^(r(-i)) 

We now return to the computation of A3 2(r, e). By redefining y^y — x]Z^z — 
yl3 2(r, e) is transformed into the following integral 

yl3 2(r,e) = 87r(n-2)A^x (A.8 

y<r ^ 

The first integral equals —^-57-- The second integral is more complicated. We will show 
in appendix D how to calculate it. In fact, it equals (cf. ( p.lO|) ) Vott j.2)_'ii ■ We obtain : 

A3,2(r,e) = 8vr2(n-2)(^^^ (A.9) 

In fact, in appendix D, the z integration is performed over the whole complex plane, while 
in the integral ( |A.8| ), the domain of integration is restricted to the disk \z\ < This 
introduces a new singularity at infinity, which we need to subtract. However, we will show 
that this singularity is cancelled by another singularity at the origin, in the computation 
that we perform in appendix D. Thus, we need only to compute the singularity at the 
origin, and this singularity must be added to the result. That the singularity at the origin 
cancels the one at infinity is obvious if we isolate these two parts in 

< e{Q)e{l)e{z)£{oo) >< £{l)e{z) > (fz 

The singularity in 2; — is / d"^ zlz]'"^'^^ = 27r-^^^^^3^ — and the one in z — > cxd is 

\z\<r/\y\ 

\z\>r/\y 
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u 



J d'^z\z\ 2 3e _ 2 -/;- (''"/ 1^1) — ^ confirming that these contributions cancel each other. 



We still need to compute the extra contribution from z 0. Some easy manipulations 
show that this contribution is the same as the one from the integral in equation ( |A.3| ). 
So the third second order term is : 



-6e~ 



A3,3(r,e) = 167r^(n-2)(^-^j (A.IO) 
Finally, in the case a = c = e and b = d = f, a straightforward calculation of the 



integral ( [A .21 ) gives us a result of the form / d'^z\z\~'^~^'^ . This does not contain any 

|z|<r 

singularities. Collecting all of our results together, we obtain the re-normalized coupling 
constant : 

9{r) = go + glA2{r,e)+g^oiAs,iir,e)+As,2ir,e)+As,3ir,e)) 

= 9o(l-9oMn-2f-^ + gl8n\n-2f-^(l + ^^^)) (A.ll) 



B Renormalisation of m 

This appendix is devoted to the computation of the renormalisation of the coupling 
constant m associated to the energy operator. As in appendix A, we will compute 
perturbatively by contracting e operators. 



B.l First order 

The first order correction to m is produced by the product of an go{ee) operator with an 
mQ^e) operator : 

\x-y\<r " 

When 6 = c, we can contract two e operators : 

gomoBi{r,e) / Yea{y)d^y 

with 

5i(r,e) = 2(n- 1) j < e{x)e{y) >q d^y 

\y—x\<T 
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= - 1) / (^) = -Mn - 1) (^) (B.l) 
B.2 Second order 

At the second order, corresponding to the product of two go{££) operators with an mo{e) 
one, we have to compute : 

J Y.^<^^^yb{x)Y.ec{y)ed{y)Y.^e{z)d^x(fy(fz (B.2) 



9l 



x-y\,\x~z\<r "^"^ 



This expression is, in fact, very similar to the one in eq.( |A.2|) . Again there will be two 
contributions. The first one corresponds to b = c d = e. We will denote this case by : 



B2,i{f,^) corresponds to y43 i(r, e) of appendix A, where the combinatorial term [n — 
2){n — 3) has been replaced by {n — l){n — 2). From eq. (|A.4|) we can read the result : 



52,i(r,e) = IQrr'in - l){n " 2) ( ^ ) (B.3) 



The second case is produced when h = c = e^a = d\ 

gomoB2,2ir, e) / Y,^a{y)d^y 



Again, -82,2 (^^j e) corresponds to ^3,2 (^^j e) computed in appendix A where {n~2) has been 
replaced by . From eq. ( [A.9| ) , we have 

-6e\ 
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S2,2(r,6)=47r2(n-1)K— 1 (B-4) 



Again, the -62,2 term must be completed by a i?2,3(?", e) term (the equivalent of 
^2,3 ('"j e)) which is : 



52,3(r,e) = 87r^(n-l)(^-^j (B.5) 

Collecting all of these results together, we then obtain for m(r) = mo(l + Bi{r,e)go + 
g'^{B2^i{r, e) + -62,2 (^^5 e) + -82,3 (?", e)) at the second order in g{r) : 

— + An\n-l)gf-:^il + ^^^ 
3e ^ 3e ^ 3e 



m(r) = mo(l - 4n{n - l)go—- + 4n\n - l)g'o^{l + —-^)) (B.6) 
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C Renormalisation of h 



We now turn to the computation of the renormahsation of the couphng constant h{r) 
associated to the a operator. Here the situation is a httle bit different, because the 
e operator and the a operator are involved. At the first order, we must compute the 
product of Qoi^s) on hQ{a). This product will give in fact no contribution (the OPA of a 
and e does not contain a operators at the first order) and so we go directly to the second 
order. 



C.l Second order 



Here, we will compute the product of two go{££) with an ho{a) operator : 

2 

Y^o J ^ea{x)eb{x)^ec{y)ed{y)^cFe{z)(fx(fy(fz 

When a = c = e^h = d the product aee will contain a a operator. By projecting aee 
on a cr(oo) operator, we obtain : 



» n 

hoglC2{r,e) I Y,(^a{z)d 

a=l 



where 



C2{r, e) = 2(n - 1) J < a{x)e{y)e{z)cr{oc) >q< e{y)e{z) >o d^yd'^z (C.l) 

\y—x\,\z—x\<r 

Using the Coulomb gas representation, this integral becomes : 

A{n-l)N. I [-^)x (C.2) 

y<r ^ ^ 

where the normalisation is the same as the one computed in the appendix A, see 



eq. (|A.7|) . The first part of this expression (the y integral) gives (— ^^), 
one is computed in appendix D, see ( p.l3| ). We obtain ; 



-6e 



C2(r,e) = (n-1)7T' 



^^P(-l)P(-i) 



while the second 



(C.3) 
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C.2 Third order 



At the third order, corresponding to the product of three goi^e) operators with a /io(c") 
one, we have to compute : 



3! 



ea{x)eb{x) ec{y)ed{y) Y £e{x)ef{x) ^ ag{z)(fx(fy(fzd/ 



u 



\x—y\,\x—z\,\x—u\<r 

Here contributions will come from the following contractions a = c = g; b = e; d = f and 
produce : 



/ n 
a=l 



where 



C3(r,e) = 4(n-l)(n-2)x (C.4) 
< a (x) e{y)e{z) a (oo) >o< e{y)e{w) >o< e{w)e{z) >o cPycfzcfw 



\y~x\,\z—x\,\w—x\<r 

Again, this expression can be computed in the Coulomb gas representation 

^ dy \ 



{n-l){n-2)N7T / \w\-^-'^'\w - ll'^-'^'d'w 



y<r 



l+9e 



U|4"i2"fc:7-T| 



where we also factorize the y and w integrations. The first two integrations will then 
produce ^'^^^^2^' while the results of the z and u integrations are given by eq. (|D.15| ) of 
appendix D. So : 



C3(r, e) = -12(n - - 2)7i^ 



9e 



l + ^(2-A) ^ ^ 
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(C.6) 



Putting together the Cj, we obtain for h{r), up to the third order : 



-6e 



hir) = h,il + in-l)gy 



3„3 



(C.7) 



-9e' 



-12{n-l)in-2)g-^n 



9e 



8 _ r^(-|)r^(i) 
,^ + 9^2-^V^H)PH) 
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D Computation of integrals 



In this appendix, we briefly present the general method to compute integrals of the 
following form 

1 = 1 \x\'-\x - - l\^''\x - y\^3dVy (D.l) 

where the integration is performed over the whole complex plane. Using the techniques 
of ref.[0] (see also this integral can be decomposed into holomorphic and antiholo- 
morphic parts : 

/ = s{b)s{b') [j+J^ + J+J^] + s{b)s{b' + 2g)J+J^ + s{b + 2g)s{b')J+J{ (D.2) 
where s{x) corresponds to sin{nx) and 

J+ = J (a, b, a', 6', g) ; = J {b, a, b\ a', ^f) 
J- = j(^b,-2- a-b-2g,b',-2- a' - b' -2g,g) 



(D.3) 



Jo 



J(-2 -a-b-2g,b,-2-a - b' - 2g, b', g) 



Here, we used the notation 

1 1 

J(a, 6, a', 6', g) = J du J dv m<^+<^'+29+1(i _ ^) V' (1 - vf3(^l - uvf (D.4) 



T{2 + a + a' + 2g)T{l + 6)1(1 + a')T{l + 2g) 



E 



-b')k{2 + a + a' + 2g)k{l + a')k 



r(3 + a + a' + 6 + +2g)T{2 + a' + 2g) fc!(3 + a + a' + 6 + +2^)^(2 + a' + 2g)k 



and 



{a)k = a{a + l)...(a + k — 1) 

The J integrals appearing in ( p.2|) are not all independent. Using contour deformation 
of integrals it can be shown that we have the following relations: 



s{2g + a + b)J{ + s{a + 6) 



s{2g + a' + b')J^ + s{a' + b')J^ 



s[a] 



s{2g + a' + b' 



(s{a')J^ + s{2g + a')4) (D.5) 



s[a 



[s{a)J:^ + s{2g + a)Jt) (D.6) 



si2g + a + b) 

With these formulas, the integrals appearing in (|A.8|) , (|C.2|) , (|C.5|) can be calculated as 
a power series of e. We will now compute explicitly these three cases : 
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We first compute the integral ( |A.8|) . The coefficients corresponding to that case 



are 



a = -2a' = 2 + 3e ; h' = -^-~e ; h=-2g = -'^-'- (D.7) 
Substituting these values in (|D.2|) , we obtain at the first order in e 



Now writing the as functions of the , with the help of relations ( p.5| , p.6| ), 
this simplifies to 

I = 3V3neJ+J+ - (D.9) 

The and J2 can now be computed explicitly using formulas (|D.3| , p.4j) . We 
obtain : 

Collecting all of these pieces together, we obtain, at the lowest order in e : 



.£1-1) 



I = VSn ^^) 3 (D.IO) 



We now compute the integral (|C2|). Here, we have : 



a = -i + Ae = -ia' ; b = 2g = b' - 1 = - e (D.U) 
Then, relating the to the J^, we obtain for / : 

3 1 

32 7re / ^_\2 32 (2 - A)7re / ^„ ^_\2 /t^ 

J=^(J2) + \ ' [J2 -Ji) +0{e^) (D.12) 



These Jj can again be computed with the help of ( |D.3| , p.4| ) : 

J J- I n-m\) , r(-|)r(i) 

^ " ^ ^ rH) -2r(-|) + r(-i) +^^^^ 



Thus : 
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The last case, corresponding to ( |U.5|) , use the following parameters 

a = -^ + Ae = -^a' ; b = 2g = - e ; b' = -^-^e (D.14) 

Performing the same manipulations, we immediately obtain 
_ 9V3ne T\-l) 3n2-A)vrer^(-ir(|) 

^ - 16 r2(-|) + 2 r2(-i) + > ^^-^^^ 
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